Bounds on isocurvature perturbations from CMB and LSS data 
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We obtain very stringent bounds on the possible cold dark matter, baryon and neutrino isocur- 
vature contributions to the primordial fluctuations in the Universe, using recent cosmic microwave 
background and large scale structure data. In particular, we include the measured temperature and 
polarization power spectra from WMAP and ACBAR, as well as the matter power spectrum from 
the 2dF galaxy redshift survey. Neglecting the possible effects of spatial curvature, tensor pertur- 
bations and reionization, we perform a Bayesian likelihood analysis with nine free parameters, and 
find that the amplitude of the isocurvature component cannot be larger than about 31% for the cold 
dark matter mode, 91% for the baryon mode, 76% for the neutrino density mode, and 60% for the 
neutrino velocity mode, at 2-a, for uncorrected models. On the other hand, for correlated adiabatic 
and isocurvature components, the fraction could be slightly larger. However, the cross-correlation 
coefficient is strongly constrained, and maximally correlated/anticorrelated models are disfavored. 
This puts strong bounds on the curvaton model, independently of the bounds on non-Gaussianity. 

PACS numbers: 98.80.Cq 



Introduction. Thanks to the tremendous develop- 
ments in observational cosmology during the last few 
years, it is possible to speak today of a Standard Model 
of Cosmology, whose parameters are known within sys- 
tematic errors of just a few percent. Moreover, the re- 
cent measurements of both temperature and polarization 
anisotropics in the cosmic microwave background (CMB) 
has opened the possibility to test not only the basic 
paradigm for the origin of structure, namely inflation, but 
also the precise nature of the primordial fluctuations that 
gave rise to the CMB anisotropies and the density per- 
turbations responsible for the large scale structure (LSS) 
of the Universe. 

The simplest realizations of the inflationary paradigm 
predict an approximately scale invariant spectrum of adi- 
abatic and Gaussian curvature fluctuations, whose am- 
plitude remains constant outside the horizon, and there- 
fore allows cosmologists to probe the physics of inflation 
through observations of the CMB anisotropies and the 
LSS matter distribution. However, this is by no means 
the only possibility. Multiple-field inflation predicts that, 
together with the adiabatic component, there should also 
be an entropy or isocurvature perturbation y|,y,|j, as- 
sociated with fluctuations in number density between 
different components of the plasma before decoupling, 
with a possible statistical correlation between the adia- 
batic and isocurvature modes |^. Baryon and cold dark 
matter (CDM) isocurvature perturbations were proposed 
long ago PI as an alternative to adiabatic perturbations. 
Recently, two other modes, neutrino isocurvature den- 
sity and velocity perturbations, have been added to the 
list [||. Moreover, it is well known that entropy per- 
turbations seed curvature perturbations outside the hori- 
zon 0, Q , so it is possible that a significant component of 
the observed adiabatic mode could be strongly correlated 
with an isocurvature mode. Such models are generically 



called curvaton models 0,0, and are now widely studied 
as an alternative to the standard paradigm. Further- 
more, isocurvature modes typically induce non-Gaussian 
signatures in the spectrum of primordial perturbations. 

In this Letter we present very stringent bounds on the 
various isocurvature components, coming from the tem- 
perature power spectrum and temperature-polarization 
cross-correlation recently measured by the WMAP satel- 
lite [9j; from the small-scale temperature anisotropy 
probed by ACBAR 0; and from the matter power spec- 
trum measured by the 2-degree-Field Galaxy Redshift 
Survey (2dFGRS) JLjJ. We do not use the data from 
Lyman-a forests, since they are based on non-linear simu- 
lations carried under the assumption of adiabaticity. We 
will not assume any specific model of inflation, or any 
particular mechanism to generate the perturbations (late 
decays, phase transitions, cosmic defects, etc.), and thus 
will allow all five modes - adiabatic (AD), baryon isocur- 
vature (BI), CDM isocurvature (CDI), neutrino isocur- 
vature density (NID) and neutrino isocurvature velocity 
(NIV) - to be correlated (or not) among each other, and 
to have arbitrary tilts. However, we will only consider 
the mixing of the adiabatic mode and one of the isocur- 
vature modes at a time. This choice has the advantage of 
restricting the number of free parameters, and takes into 
account the fact that most of the proposed mechanisms 
for the generation of isocurvature perturbations lead to 
only one mode. The first bounds on isocurvature pertur- 
bations assumed uncorrelated modes [13. but recently 
also correlated ones were considered [I3lll4 ITEl Il6| . 

The present analysis neglects the possible effects of 
spatial curvature, tensor perturbations and reionization. 
Therefore, each model is described by nine parameters: 
the cosmological constant Q.\, the baryon density cjb = 
fie^ 2 , the cold dark matter density w c d m = ^cdm/i , the 
overall normalisation A, the isocurvature mode relative 
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amplitude a and correlation f3, the adiabatic and isocur- 
vature tilts (n a( j, n.; so ), and finally a free bias 6 associated 
to the 2dF power spectrum. We generate a 5-dimensional 
grid of models {A, a, [3, and b are not discretized) and 
perform of Bayesian analysis in the full 9-dimcnsional 
parameter space. At each grid point, we store some Ci 
values in the range < I < 1800 and some P(k) values in 
the range probed by the 2dF data. The likelihood of each 
model is then computed using the software or the detailed 
information provided on the experimental websites, using 
1398 points from WMAP, 11 points from ACBAR and 32 
points from the 2dFGRS. For parameter values which do 
not coincide with grid points, our code first performs a 
cubic interpolation of each power spectrum, accurate to 
better than one percent, and then computes the likeli- 
hood of the corresponding model. From the limitation of 
our grid, we impose the flat prior ri; so > 0.6 (as preferred 
by inflation). The other grid ranges are wide enough in 
order not to affect our results. 

For the theoretical analysis, we will use the notation 
and some of the approximations of Ref. ^4| . For instance, 
the power spectra of adiabatic and isocurvature perturba- 
tions, as well as their cross-correlation, are parametrized 
with three amplitudes and two spectral indices, 
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where kg is an arbitrary pivot scale. We also assume 
that the correlation coefficient cos A is scale-independent 
(Note that in Ref. ^(| this assumption was relaxed). 
In order to evaluate the temperature and polarization 
anisotropies, one has to calculate the radiation trans- 
fer functions for adiabatic and isocurvature perturbations 
and compute the total angular power spectrum as [T^j 

Ci = Cf + B 2 Cf ° + 2B cos A C, corr , (2) 

where B is the entropy to curvature perturbation ratio 
during radiation domination, B = S/1Z. We will use 
here a sli ghtl y different notation, used before by other 
groups 0,ll3|, where 

Q = (1 - a) Cf d + a Cf ° + 2/3 y/a(l-a) Cf°" . (3) 

The two notations are related by 

a = B 2 /(l + B 2 ) , /3 = cosA. (4) 

This notation has the advantage that the full parameter 
space of (a, 2/3-%/ a(l — a)) is contained within a circle 
of radius 1/2. The North and South rims correspond 
to fully correlated (f3 = +1) and fully anticorrelated 
((3 = — 1) perturbations, with the equator correspond- 
ing to uncorrelated perturbations (/? = 0). The East 



and West correspond to purely isocurvature and purely 
adiabatic perturbations, respectively. Any other point 
within the circle is a particular admixture of adiabatic 
and isocurvature modes. 

In order to compute the Ci coefficients of the tempera- 
ture and polarization power spectra, as well as the matter 
spectra P(k), we have used a CMB code developed by one 
of us (A.R.) which coincides, within 2% errors, with the 
values provided by CMBFAST for AD and CDI modes, 
and also includes the neutrino isocurvature modes as well 
as the cross-correlated power spectra. Note that the code 
defines the tilt of each power spectrum with respect to 
a pivot scale fco corresponding the present value of the 
Hubble radius, while CMBFAST uses k = 0.05 Mpc" 1 : 
so, the comparison of our results with those of Refs. 
for the CDI mixed model is not straightforward. 

TABLE I: The one-dimensional 2-a ranges on the isocurvature 
mode coefficients for the various models, uncorrelated (middle 
column) and correlated (right column). 
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CDI 


< 0.31 


< 0.47 


-0.31 


to 


0.31 


BI 


< 0.91 


< 0.95 


-0.80 


to 


1.00 


NID 


< 0.76 


< 0.77 


-0.30 


to 


0.78 


NIV 


< 0.60 


< 0.60 


-0.77 


to 


0.35 



Results. We now describe the different bounds, sum- 
marized in Table I. For each specific mode, the probabil- 
ity distribution for a (uncorrelated case) and the likeli- 
hood contours in the (a, (3) plane (correlated case) are 
shown on Fig. 1. The best-fitting parameter set for each 
model is given in Table II. 

CDM isocurvature. Assuming no correlation between 
the adiabatic and CDI modes, we find an upper bound 
a < 0.31 at 95% c.L. The deviation from a pure adia- 
batic model does not improve the goodness-of-fit, since 
the minimum \ 2 value goes from 1478.8 to 1478.3, while 
the number of degrees of freedom decreases from 1435 to 
1433. The situation does not change significantly when 
we include a possible correlation: the minimum \ 2 is still 
around 1478.0, and the individual 2-er bo unds on th e coef- 
ficients read a < 0.47 and -0.31 < 2/3 v /a(l - a) < 0.31 
at 95% c.L. Most of the constraints on this model come 
from the WMAP TT power spectrum. For comparison, 
we repeated our analysis replacing WMAP and ACBAR 
by the Wang et al. [17| CMB data compilation, which 
includes all the pre- WMAP temperature data. In that 
case, the results exhibit a wide parameter degeneracy (as 
indicated by previous studies [l3l H3 | ) , and a large frac- 
tion of anti-correlated isocurvature modes cannot be ex- 
cluded, up to (a, (3) = (0.9, —1), provided one allows for 
a very large baryon fraction. Therefore, we conclude that 
the WMAP TT data is very powerful in constraining the 
isocurvature fraction. On the other hand, we checked ex- 
plicitly that the effect of the WMAP TE power spectrum 
is not very strong with respect to that of TT. 

In the curvaton scenario, the case in which the CDM 
is created before the curvaton decay while the curva- 
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FIG. 1: a) the likelihood function of the isocurvature fraction a, for three different types of uncorrelated isocurvature modes 
(i.e., with the prior f3 = 0); b) the 1 and 2-a contours of a and the cross-correlated mode coefficient 2(3^/ q(1 — a), for the 
CDM isocurvature mode: the small (red) contours are based on all the data, with one flat prior ni so > 0.6, while the large 
(green) ones show the situation before WMAP, with an additional prior lob < 0.037; c) same as b) for NID; d) same as b) for 
NIV. The marginalization is approximated by a maximum likelihood fit of the other parameters for each pair of values (a, f3). 



ture perturbations are small, corresponding to (a, (3) = 
(0.9, —1) in our notation Q), is completely excluded, as 
already emphasized in Ref. [15|. The case in which the 
CDM is created by the decay of the curvaton leads to 
f3 = 1; with such a prior, we obtain a sharp 2-a bound: 
a < 0.04 or, in the notation of Ref. 14], B < 0.2 
(the authors of Ref. |15| obtain a more restrictive bound 
B < 0.43 ^b/^m ~ 0.1, presumably because they use 
the curvaton-motivated assumption ni SO = n ac i). We 
can rewrite this bound as a constrain on the fraction 
r = il (Ti decay since, in these models, B = 3(1 — r)/r. In 
our case, B < 0.2 implies 1 — r < 0.0625. 

Baryon isocurvature. The case of baryon isocurva- 
ture modes is qualitatively similar to that of CDI modes, 
since the spectra are simply rescaled by a factor ^i/^cdm 
(SlB/^cdm) for the isocurvature power (cross-correlation) 
components: thus, significantly larger values of a will be 
allowed in the baryon case. Some approximate results 
for the BI modes could be deduced from the CDI results 



TABLE II: The best fit values of the parameters for the adia- 
batic (AD), uncorrelated (CDI, etc.), and correlated (c-CDI, 
etc.) isocurvature models, with the corresponding % 2 an d 
number of degrees of freedom v. For the NIV uncorrelated 
model, the best-fit occurs for a — (i.e., purely adiabatic). 



model 


wb Wcdm n a d rii so a P 


X> 


AD 


0.021 0.12 0.70 0.95 - 


1478.8/1435 


CDI 
BI 
NID 

NIV 


0.023 0.12 0.73 0.99 1.02 0.10 - 
0.023 0.12 0.73 0.99 1.02 0.72 - 
0.023 0.12 0.73 0.99 0.95 0.37 - 
0.021 0.12 0.70 0.95 - - 


1478.3/1433 
1478.3/1433 
1478.2/1433 
1478.8/1433 


c-CDI 
c-BI 
c-NID 
c-NIV 


0.022 0.12 0.71 0.97 1.23 0.001 1 
0.022 0.12 0.71 0.97 1.23 0.03 1 
0.022 0.12 0.73 0.97 1.04 0.10 0.26 
0.021 0.12 0.71 0.95 0.71 0.03 -1 


1478.0/1432 
1478.0/1432 
1477.7/1432 
1477.5/1432 



by an overall rescaling of the a parameter. However, 
we performed an exact analysis, and found the following 
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bounds: a < 0.91 (uncorrelated case), or a < 0.95 and 
—0.8 < 2(3^Ja{\ — a) < 1 (correlated case). In each case, 
the minimum \ 2 is, by definition, the same as for the CDI 
case. We provide the best-fit parameters in Table I, but 
for brevity we do not show this case in the figures. 

In the curvaton scenario, the case in which the baryon 
number is created before the curvaton decay while 
the curvature perturbations are small, corresponding to 
(a, j3) — (0.9, —1), is also excluded at several cr's, as 
previously found by in Ref. [l5j |. The case in which the 
baryon number is created by the decay of the curvaton 
predicts (3 = 1; we then find a < 0.5, or B < 1. 

Neutrino density isocurvature. The 2-cr bound for the 
uncorrelated case is a < 0.76, while for the correlated 
adiabatic and neutrino mode we obtain a < 0.77 and 
-0.30 < 2/3y / a(l - a) < 0.78. The best-fitting models 
for the two cases are still extremely close to the adia- 
batic model, leading to no significant improvement in 
the likelihood, A% 2 = 0.6 vs. A% 2 = 1.1, respectively. 
Note that with the pre-WMAP CMB data plus the 2dF- 
GRS power spectrum, significantly larger fractions, up to 
(a, (3) = (0.9, —0.6), of correlated NID modes were then 
allowed, but are now excluded with great confidence. 

The authors of Ref. @ discuss an interesting mecha- 
nism under which the curvaton could generate some fully 
correlated or anticorrelated NID modes, through inho- 



mogeneities in the neutrino/antineutrino asymmetry pa- 
rameter. However, our bounds for these cases are very 
restrictive, a < 0.18 (correlated), or a < 0.02 (anticorre- 
lated) . In the notation of Ref. [lj] , this is equivalent to 
B < 0.47 or B < 0.14, respectively. 

Neutrino velocity isocurvature. Again, the presence of 
such a mode is not significantly favored: in the uncorre- 
lated case, the best-fitting model is still the pure adia- 
batic one (a = 0), while in the correlated case the \ 2 i m_ 
provement is only A% 2 = 1.3. The impact of WMAP on 
these bounds is spectacular. Using only the pre-WMAP 
CMB data compilation and the 2dFGRS power spec- 
trum, we find that half of the parameter space (with 
< a < 1 and < (3 < 1) was allowed until recently 
(at the expense of an unnaturally high baryon fraction 
and scalar tilt). 

Conclusions. Using the recent measurements of tem- 
perature and polarization anisotropies in the CMB by 
WMAP, as well as the matter power spectrum measured 
by 2dFGRS, we obtained very stringent bounds on a pos- 
sible isocurvature component in the primordial spectrum 
of density and velocity fluctuations. We have considered 
both correlated and uncorrelated adiabatic and isocurva- 
ture modes, and put strong constraints on the curvaton 
scenario of primordial perturbations. 
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